Single-Cycle Gap Soliton in a Subwavelength Structure 
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We demonstrate that a single sub-cycle optical pulse can be generated when a pulse with a few 
optical cycles penetrates through resonant two-level dense media with a subwavelength structure. 
The single-cycle gap soliton phenomenon in the full Maxwell-Bloch equations without the frame of 
slowly varying envelope and rotating wave approximations is observed. Our study shows that the 
subwavelength structure can be used to suppress the frequency shift caused by intrapulse four-wave 
mixing in continuous media and supports the formation of single-cycle gap solitons even in the 
case when the structure period breaks the Bragg condition. This suggests a way toward shortening 
high-intensity laser fields to few- and even single-cycle pulse durations. 

PACS numbers: 42.65.Re, 42.65.Tg, 42.50.Gy 



The recent development of ultrafast science technol- 
ogy allowed the generation of light pulses with durations 
down to few optical cycles. There is an increasing number 
of investigations and experiments which involve such ul- 
trashort pulses in femtosecond and attosecond domains 
Applications also include control of chemical reac- 
tions in physicochemical processes [2j, femtochemistry 
0| , biological multiphoton imaging [j] and coherent con- 
trol schemes In order to study light-matter interac- 
tion under extreme conditions as the pulse approaches 
the duration of a single optical cycle, great efforts for 
the generation of still shorter pulses have been made 
[§]. In particular, for continuous media the concept of 
single-cycle nonlinear optics @ as well as the notion of 
extreme nonlinear optics (Io| were introduced. For these 
processes, the traditional frame of the slowly varying en- 
velope approximation (SVEA) and the rotating wave ap- 
proximation (RWA) is invalid [111 ]. 

In the last decade, there has been a rapid progress in 
the field of materials with periodic structures, such as 
Bragg grating, photonic band gap crystals and waveg- 
uide arrays [121 ] . Such systems exhibit qualitatively 
novel and fascinating linear-optical, nonlinear-optical 
and quantum-optical properties which provides an at- 
tractive way to control the light propagation and light- 
matter interaction. Nonlinear pulse propagation in pe- 
riodic structures is of both fundamental and applied in- 
terest, particularly, as a potential basis for nonlinear fil- 
tering, switching and distributed-feedback amplification 
[l3j |. Few-cycle pulse propagation in continuous media 
can yield solitons. Solitons can also exist in periodic 
structures which is usually referred to as gap soliton since 
they only exist within the forbidden gap 14|,[l5|. The gap 
soliton was primarily found in periodic Kerr-nonlinear 
media. An essentially different mechanism of gap soli- 
ton generation with self-induced transparency has been 
revealed in Bragg grating which consists of a periodic 
array of thin layers of a resonant two-level medium sep- 
arated by half- wavelength nonabsorbing dielectric layers 
TH 17j |. It should be pointed out that for such a struc- 



ture, the thickness of each two-level medium layer is much 
smaller than the wavelength of the electromagnetic wave 
which propagates. Optical soliton propagation results 
from an interplay between dispersion and the nonlinear 
response of the medium which balance the nonlinear and 
linear effects [HI. Optical temporal solitons have be- 
come a promising candidate for optical communication 
networks. Thus, it is an important issue whether or not 
optical solitons can exist under extreme conditions, such 
as pulses with few optical cycles or less formed in a pe- 
riodic structure. The ability to reduce the cycle number 
in these systems would be of particular significance. 

Therefore, in this Letter we present theoretical evi- 
dence of a sub-cycle gap soliton. An extremely short 
pulse propagating through a one-dimensional periodic 
structure that consists of periodic thin layers contain- 
ing resonant two-level atoms separated by a transparent 
material is considered (see Fig.fTJ). We assumed that the 
extremely short pulse propagates along the £ axis to an 
input interface of the periodic structure at £ = (a set 
of dimensionless time and space coordinates were taken, 
namely T = tujQ and £ = zj Ao with the carrier frequency 
wo and wavelength Ao). First, the incident pulse moves 
in the free space and then partially penetrates into the 
subwavelength structure. In the following, the penetrat- 
ing part passes through and exits into the free space. The 
length of the structure is L. We found that at resonance 




FIG. 1: (color online) Spatial distribution of atomic medium. 
The red curve illustrates the incident wave with few optical 
cycles. Symbol "— s>" denotes the propagation direction axis. 
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and for a periodic structure satisfying the Bragg condi- 
tion for moderate densities, the penetrating and reflected 
pulse spectrum is located around the transition frequency 
which is distinct from the results of few-cycle propaga- 
tion in a continuous media. There, a large frequency shift 
of the transmitted and reflected pulses was reported [1] . 
For particular parameters, we obtain a gap soliton with 
less than one optical cycle even when the structure period 
does not fulfill the Bragg condition. Thus, the system we 
propose allows the observation of different tendencies in 
the spectral and temporal evolution of ultrashort pulses 
leading to intresting and practically significant regimes 
of ultrashort-pulse quantum dynamics. Finally, these 
periodic structures can be realized by using a periodic 
doped fiber waveguide [IH or quantum wells embedded 
in a semiconductor structure |19| . 

In order to investigate the electromagnetic field prop- 
agation in the periodic structure, we adopt the full 
Maxwell-Bloch equations solved by a finite difference 
time-domain method [20j | which has been proved to be an 
accurate tool. In this case, the evolution of the few-cycle 
electromagnetic field amplitudes E x (T, £) and H y (T, £) , 
the complex atomic polarization ^u(T, £) + |i>(T, £), and 
the inversion population u>(T, £) obey the following set 
of Maxwell-Bloch equations: 
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Here, E = E x /Eq and H = y/ fi / 'eoH y / 'Eq are the nor- 
malized dimensionless components of the few-cycle pulse 
(Eq corresponds to the peak amplitude of the electric 
component of the incident pulse) . fin = gEo / h with the 
dipole moment g is the peak Rabi frequency of the in- 
cident pulse. T\ and T 2 are, respectively, the lifetime of 
the excited state and the dephasing time. In the above 
equations, we have assumed that the transition frequency 
of the two-level atom is equal to the carrier frequency u>q 
of the pulse. Further, w c (£) = N(£)g 2 /e h with N(£) 
denoting the density distribution of the atomic medium 
and which is given by the expression 




if £ e [2nS, 2nS + 5), 

if £ £ [2nS + S, 2nd + 26), 



(2) 



where d and 8 express the atom density and the thickness 
of one layer, n is an integer between [0, LJ (28) — 1] . The 
density distribution of the atomic medium is shown in 
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FIG. 2: (color online) Electric field profiles for different time 
moments, (a) t = 0.1 ps; (b) t — 0.3 ps; (c) t — 0.5 ps; (d) 
t — 0.7 ps. The other atomic and structure parameters are: 
fi = 1.4 fiT 1 , d = 0.2, 5 = 0.25 and £ £ [0,200]. (a) and 
(d) correspond to the electric field profiles near the input face 
and outside the exit face (vertical line in d), respectively. 



Fig.HJ The incident pulse we consider is a few-cycle pulse 
with a hyperbolic secant envelope, i.e., il(T — 0, £) = 
fi cos[27r(£— £ )]sech[1.76x27r(£— £ )/Tp]) where r p is the 
full width half maximum (FWHM) of the pulse intensity 
envelope (see the red curve in Fig. [T]) . The choice of £ 
ensures that the pulse locates in the free space and the 
partial pulse in the structure can be neglected at t = 0. 
The two-level medium is initialized with u = v — and 
w = —1 at t = 0. The material and laser parameters 
are chosen as follows: r p = 5 fs, u>o = 2.3 fs _1 (An = 
830 mn), g = 2 x 1CT 29 A s m, Ti = 1 ps, f 2 = 0.5 ps, 
and N = 4.4 x 10 20 cm~ 3 . The above parameters give 
an atomic density d = 0.2. The peak Rabi frequency 
fin = 1.4 fs -1 corresponds to an input envelope area 
A = f2nT p 7r/1.76 = 47T. The cycle number of the initial 
pulse is about 1.83. 

We have obtained numerical results of Eqs. (fTJ) for dif- 
ferent periodic structures. Firstly, let us investigate the 
few-cycle pulse propagation in a resonant structure in 
which the structure period 28 exactly satisfies the Bragg 
condition: 28 = m/2 (the dimensionless coordinate is 
taken and m is an integer here). In our case, the period of 
the structure is set to be equal to 1/2. The electric field 
profiles of the penetrating pulse inside and outside the 
periodic structure for different time moments are shown 
in Fig. [2j As can be observed, the pulse can maintain the 
same envelope to propagate which is a typical behavior of 
the gap solitons in periodic structures. Remarkably, the 
cycle number of the transmitted pulse is reduced in com- 
parison to the initial pulse. We obtain a 0.85 cycle for the 
output pulse in Fig.[2][d). To further describe the system, 
we plot in Fig. |3] the inversion as well as the electric field 
profile at particular time-space coordinates. The distri- 
bution of inversion w is not continuous and in certain 
layers one can create population inversion. Therefore, 
the sideband parts of the propagating pulse annihilate 
while the central part slightly amplifies. Evidently, the 
output pulse energy is not increased in comparison to 
that of the initial applied pulse. 
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FIG. 3: (color online) Inversion population distribution (red 
curve) and electric field profile (blue curve) at time t = O.lps. 
flo = 1.4 fir 1 , 5 = 0.25, d = 0.2 and £ £ [0, 108]. 



Subsequently, in order to clarify the soliton propaga- 
tion in detail, we present in Fig. @] the evolution of the 
pulse spectrum during propagation through a periodic 
structure when 5 = 0.25. Inspecting Fig. [4j we observe 
that the reflected and the penetrating spectra are sit- 
uated around the input central frequency ujq which is 
distinct from the results of few-cycle propagation in the 
continuous media ||. There, for d = 0.2 and £ = 108, the 
carrier frequency of the output pulse locates at 1.16wo- 
In our case, however, the large frequency shift of the 
transmitted and reflected pulses caused by intropulsc 
four-wave mixing can be efficiently suppressed. Peri- 
odic structures can usually offer optical band gaps simi- 
lar to the flow of electrons in semiconductors and modu- 
late the transmission of waves with different frequencies. 
In the same way, the periodic structure in our model is 
able to modify the dispersion property of the two-level 
atomic medium and destroy the phase-matching condi- 
tion of four-wave mixing process which guarantees that 
the deformation of the penetrating pulse does not hap- 
pen in our model. From here, we can conjecture that the 
periodic structure provides an alternative method for ac- 
curately preparing few-cycle pulses which will allow one 
to approach the experimental study of single-cycle optics 
[§]. Thus, we have demonstrated that single sub-cycle 
gap solitons can be formed in a periodic structure. 

In the following, we will discuss the case of a few-cycle 
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FIG. 4: (color online) Evolution of the pulse spectrum during 
propagation inside and outside the periodic structure. The 
atomic and structure parameters are the same as in Fig. ©. 
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FIG. 5: (color online) Electric field profiles for different time 
moments and structure periods near the input face [(a),(b)] 
and outside the exit face [(c), (d)]. (a) t = O.lps, S = 0.25; 
(b) t = O.lps, 5 = 0.1; (c) t = 0.4ps, S = 0.25; (d) t = 0.4ps, 
8 = 0.1. Red line on the left of (a) and (b): incident pulse at 
t = 0. Here fl = 1.4 fir 1 , d = 0.2 and £ £ [0, 108]. 



pulse propagation in a tighter periodic structure which 
violates the Bragg condition. As a comparison, in Fig. [SJ 
the initial, reflected, penetrating and transmitted pulses 
near the input face and outside the exit face are shown 
for different structure periods. In the simulation, the pa- 
rameters are the same except for the structure period 28. 
For the case of tighter periodic structure (5 = 0.1), the 
penetrating part is strong enough to split into two pulses 
where the stronger of which moves faster than the weaker 
one that is unstable and decays quickly (see Fig. and 
Fig. [5jl). More energy of the incident pulse is located 
in the structure in this situation which is different from 
the case of S = 0.25 (compare Fig. [5^ with Fig. [SJd). 
At the same time, by comparing with the profiles of the 
penetrating and transmitted pulses, one can observe that 
the incident pulse is dramatically shortened. We obtain 
the cycle number, N c , for the FWHM of the intensity of 
the output pulse in Fig. [5jc,d) as N c ~ 1.76r/(27r) with 
r w 2.5. Thus, the cycle number iV c is about 0.7 which 
indicates that the gap soliton with less than one optical 
cycle can be formed in a subwavelength structure. The 
pulse propagation for other cases (i.e. for longer pulses 
and structure periods, slightly different and random layer 
widths, etc.) has been simulated, too, and the behavior 
of the penetrating pulse is similar to the situations shown 
in Fig. [5]- Fig. though for longer pulses one can obtain 
fewer output stable pulses. If the atomic density is in- 
creased, the single-cycle gap soliton cannot be obtained 
for larger structures. The reason for this is that the re- 
flection on the interfaces of the atomic layers is enhanced 
rapidly with the increase in the atomic density. The out- 
put field spreads and becomes weak. Similar things can 
happen if one moderately increases the initial pulse in- 
tensity. However, shorter structures (£ ~ 20) can lead to 
an output single-cycle pulse, too. If 0.2 < d < 1, its car- 
rier frequency slightly shifts towards higher frequencies 
when Oq approaches ujq (the shift is smaller than that in 
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FIG. 6: (color online) Electric field profiles and corresponding 
fitting curves in space and time coordinates (a) for t — O.lps; 
(b) for £ = 30. The black and red lines represent the nu- 
merical result and the analytic solution fl + given in Eq. ([3|, 
respectively. Q = 1.4 fir 1 , S = 0.25, d = 0.2 and £ G [0, 200]. 



a continuous sample with identical parameters). 

In general, gap solitons are generated as a result of 
interaction of counterpropagation waves. As remarked 
above, the existence of gap soliton solution of the en- 
velope equations is restricted by the validity of SVEA 
and RWA which is violated for the conditions consid- 
ered here. However, in the above two approximations, 
the pulse propagation in a resonant periodic structure 
can be described by a set of semiclassical coupled-mode 
Maxwell-Bloch equations. Then, assuming that the pulse 
consists of two strong coupled modes that propagate in 
the structure, one can obtain a sec/i-form soliton solution 
which has the following form 17J: i? ± = Ml^ /g, where 



fi±sech[(27r£ - v g T)/(v g r)] cos(2< =f T + 4>). (3) 



Here, E + and E" denote the electric-field amplitudes of 
the forward and backward waves. v g and r are group ve- 
locity and pulse duration, respectively. 4> is an arbitrary 
phase factor while ilf^ — u>o(l ± l/v g )/r. The relation- 
ship of parameters v g and t is d 9 2 = 1 + 2uj £ t 2 S/ojo 
with S depending on the broadened line form [17[ • With 
the help of these results, one obtains the envelope ar- 
eas of the forward and backward waves: A ± = J ^(1 ± 
^ 1 )sech[(2 7 r^ - v g T)/(v g T]]dt = n(l ± v' 1 ) which is re- 
ferred to as gap 27T soliton [ll [l7| . If the field amplitude 
is strong enough, the cycle number of the stable soliton 
will be small or less than one. 

As an example, in Fig. [6l the electric field profiles 
and corresponding fitting curves in space and time co- 
ordinates are illustrated. Here, the group velocity v g = 
(£2 — £i)/(^2 — ti) ~ 0.85 where £2 and £1 are the po- 
sitions of the pulse maxima at i2 and t\, respectively, 
(see, for example, Fig. [2a, c) whereas velocity unit is the 
light velocity in vacuum). The pulse duration is r w 2.5. 
Due to the fact that the pulse curves are mainly deter- 
mined by the forward wave in the early stage of the for- 
mation of gap soliton, the fitting function is taken to be 
of the form of a forward wave f2 + , i.e., Vt + — (loq/t)(1 + 
l/v ff )sech[27r(£ - £m)/(« s r)] cos (2tt(£ - £ m ) + ^) when 
t is fixed while D,+ = (w /t)(1 + l/u 9 )sech[-(T - 



T m )/r] cos (— (T - T m ) + <j> T ) when £ is fixed. {£ m , T m } 
and 4>t} describe the corresponding positions of the 
maxima of the profile pulse given by Eqs. ((T|). In partic- 
ular, we obtain £ m ps 16.6, T m w 325.5, <f>^ ~ 0.877T and 
(j>T ~ 1-97T for the curves shown in Fig. El Here, one inter- 
esting result is that the sech-iorm solution derived from 
the semiclassical coupled-mode Maxwell-Bloch equations 
within the frame of SVEA and RWA coincides well with 
the numerical results in our model. With the further evo- 
lution of the penetrating pulse, the electric field profiles 
will depend on the forward and backward waves. 

In summary, we demonstrated the formation of a single 
sub-cycle gap soliton when a few-cycle laser pulse pene- 
trates a subwavelength structure. The carrier frequency 
of the output pulse is located around the frequency of 
the initial laser pulse. The proposed scheme works even 
for the case of broken Bragg condition as well as slightly 
different random layer widths. The corresponding ana- 
lytical solution of the single-cycle gap soliton was given. 

We acknowledge helpful discussions with Christoph H. 
Keitel and Karen Z. Hatsagortsyan. 
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